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Abstract 



In this paper we study, Prob(n, a, 6), the probabihty that all the eigenvalues of 
finite n unitary ensembles lie in the interval (a, b) . This is identical to the probability 
that the largest eigenvalue is less than b and the smallest eigenvalue is greater than 
a . It is shown that a quantity allied to Prob(n, a, b) , namely, 



Hn{a,b) :-- 



d_ d_ 
da db 



In Prob(n, a, b), 



in the Gaussian Unitary Ensemble (GUE) and 



Hn{a,b) :-- 



a— + b— 
da db 



In Prob(n, a, b), 



in the Laguerre Unitary Ensemble (LUE) satisfy certain nonlinear partial differential 
equations for fixed n, interpreting Hn{a,b) as a function of a and b. These partial 
differential equations maybe considered as two variable generalizations of a Painleve 
IV and a Painleve V system, respectively. As an application of our result, we give an 
analytic proof that the extreme eigenvalues of the GUE and the LUE, when suitably 
centered and scaled, are asymptotically independent. 
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1 Introduction 



In the theory of random matrices, the study of eigenvalue distribution attracts the most in- 
terest and has many apphcations in both mathematic, physics and wireless communications; 
see for example |3l UHl |211 l2Z]- It is by now a classic result that the largest eigenvalue dis- 
tribution of the Gaussian Unitary Ensemble (GUE) and Laguerre Unitary Ensemble (LUE), 
the celebrated Tracy- Widom II distribution [21], denoted as F2 , is given by a integral in- 
volving the Hastings-McLeod solution of the Painleve II equation. For the GUE, it can be 
immediately seen via a change of variables that the analogous result holds for the smallest 
eigenvalue. These distributions emerge after centering the extreme eigenvalues at ±A/2n, 
the edges of the GUE spectrum, followed by a scaling with respect to the density at the soft 
edge; see [21] for the original derivation of this particular P77 . In [30] a finite n version of 
the result was found, which turns out to be a Piy , with the limiting F2 after the above 
centering and scaling of the largest eigenvalues. For LUE, the distribution is obtained after 
centering the largest eigenvalues at 4n and a scaling at the soft edge of the LUE spectrum. 
The finite n analogue in this case was found to be a Py [30] . 

In this paper we consider a related problem. We are concerned with the probability that all 
the eigenvalues are in an interval (a, h) . This is of course equivalent to the probability that 
the largest eigenvalue is less than h and the smallest one greater than a . It turns out that 
this probability is related to another expression that is the solution of a nonlinear partial 
differential equation (pde). In the GUE case, this pde maybe interpreted as a two variable 
version of a Painleve IV cr -form, since in the limit a fixed, 6 — )■ 00 or 6 fixed, a — )■ —00 , 
the pde reduces to the ode corresponding to the "left" or "right" edge version of Pjy . We 
have a similar interpretation for the pde derived in the LUE two variable version 

of a Painleve V a -form. See also [2j where the authors derived pdes for the logarithm of 
this probability which are related to KP equations using vertex operators and the associated 
Virasoro constraint [T]. 

Our method is based on a ladder operators formalism for orthogonal polynomials [13] and 
the associated compatibility conditions. This is, by now, a well-known method that has 
been applied to obtain exact solutions in a classical case [E] and adapted to orthogonal 
polynomials on the unit circle [1] intimately related to the theory of Toeplitz matrices. One 
may find in [TB] a list of references on this formalism. Recent examples on the applications 
of the ladder operators with the resulting Painleve equations can be found in [51 [TTl [T21 [IZJ 
[T^ [2U| [22] ; see also pK] for an application to the information theory of multiple input and 
multiple output wireless communications which involves certain deformation of the LUEs. In 
particular, a comparison between the ladder operator theory and the isomonodromy theory 
of Jimbo-Miwa-Ujimo [26] is carried out in [16] and [22] for different specific Hermitian 
random ensembles. The extension of ladder operators to discrete orthogonal polynomials 
and g -orthogonal polynomials is given in [24] and [15J, respectively. From these extensions. 
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it is shown that the recurrence coefficients of certain discrete or g -orthogonal polynomials 
are related to the discrete or g-Painleve equations, respectively. We refer to [H El El [23] for 
investigations of this aspect. 

Recently, based entirely on the estimate on the integral operator in a Fredholm expansion, 
the authors of [6] showed that the extreme eigenvalues of GUE, when suitably centered and 
scaled, are asymptotically independent random variables. As an application of our result, we 
give a proof that is an analytic counter-part of this probabilistic version. Our strategy is to 
scale the two variables in the pde obtained near both edges of the GUE spectrum, which will 
lead to a limiting pde. It turns out that the solution to the limiting pde has a particular form 
which implies that the joint probability density is a product of two independent densities, 
that is, there is asymptotic independence. With different and delicate scalings, we show 
that the solution of limiting pde in LUE case is also asymptotically equal to the sum of 
Tracy- Widom left and right distributions. 

The rest of this paper is organized as follows. In section 2, we give a summary of the 
ladder operators method, the associated compatibility conditions and a summation identity. 
From these conditions, we derive a system of non-linear difference equations satisfied by 
auxiliary variables that appear naturally in this approach. These difference equations will 
be instrumental in our further derivation of the pdes. We study GUE in section El and 
LUE in section |H As an application of our result, we give an analytic proof that the 
extreme eigenvalues of GUE and LUE, when suitably centered and scaled, are asymptotic 
independent. 



2 Preliminaries 



In the theory of Hermitian random matrices, one encounters the following (un-normalized) 
joint probability density of eigenvalues {xj}^^i : 

71 

p{Xi,...,Xn) = n [^n{x)]^Y[w{Xk), (2.1) 

l<i<j<n k=l 

where 

An(a;):= n (^^-^J-), (2.2) 

l<j<i<n 

and w is a weight function defined on an interval ICR. The GUE corresponds to 
w{x) = and / = M, while the LUE corresponds to w{x) = x°'e~^ , a > and 

/ = (0, oo) . In what follows everything depends on n , but the dependence will only be 
displayed when it is easy to do so or needed for clarification. 
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Denote by Prob(r;,, a, b) the probability that all eigenvalues lie in an interval (a, b) G I . We 
then have 

It is a well-known fact (cf. [28]) that the multiple integrals in fl2.3p can be expressed as 
determinant of Hankel matrix generated by the weight function w . The moments of the 
weight are defined by 

f' ■ 

fij{a,b) := / x^w{x)dx, j = 0, 1,..., (2.4) 



and the Hankel determinant is the determinant of the moment matrix (nj+k)^ f^Lo '■ 

Dnia, b) := det (fXj+kia, b)Y.~^^Q 

n~l 



j=0 



1 r " 

= - / [A,,ix)]^T\w{xk)dxk. (2.5) 

J{a,b)" t=i 

Here, hi{a,b) is the square of the norm of the (monic) polynomials orthogonal with 
respect to w over (a, 6) , i.e.. 



6 

Pi{x, a, b)Pj{x, a, b)w{x)dx = hi{a, b)Sij. (2.6) 



The monic polynomials Pn{z,a,b) are normalized as 

Pn{z, a, b) = z'' + pi(n, a, 6)^""^ + ... + P^O, a, b). (2.7) 

Therefore, 

Prob(n,a,6) = ^-^ = n^. (2.8) 

It is then clear that in order to compute the probability Prob(n, a, b) , we need to compute the 
product of the norms /ij(a, b) and thus we need information about the orthogonal polynomi- 
als. According to the general theory of orthogonal polynomials, an immediate consequence 
of (12.61) is the three terms recurrence relations: 

zPn{z, a, b) = Pn+i{z, a, b) + a„(a, b)Pn{z, a, b) + /3„(a, b)Pn-i{z, a, b) (2.9) 

with the initial conditions 

Po[z) = l, /3oP-i(^)=0. (2.10) 



4 



An easy consequence of the recurrence relation is 

an{a, h) = Pi(n, a, b) - p^{n + 1, a, h). (2.11) 

In the next two sections, we will give an account for a recursive algorithm for the determi- 
nation of the recurrence coefficients a„ , /3„ in the GUE and LUE, respectively, based on a 
pair of ladder operators and the associated supplementary conditions. It will become clear 
that the determination of a„ and /3„ will yield the necessary information to find conditions 
on hn{a, b) . The following three lemmas can be found in, for example, [16j and the extensive 
references therein. For convenience, we suppress the a , b dependence in a„ , /3„ and hn 
in these lemmas. 

Lemma 1. Suppose that w is such that the moments 

Hi{a,b) := / x''w{x)dx, z = 0,l,... (2.12) 



exist and that v = — In w has a derivative in some Lipshitz class with positive exponent. 
The lowering and raising operators satisfy the following: 

P'^{z, a, b) = -Bn{z, a, b)Pn{z, a, b) + /3„v4„(2;, a, 6)P„_i(z, a, 6), (2.13) 
P;;„i(z,a,6) = [Bn{z,t) + \/{z)\Pn-i{z,a,b) - An^i{z,a,b)Pn{z,a,b), (2.14) 

where 

w{y)P^{y,a,b) ' ^ 



An{z,a,b) :-- 



hn Ja Z-y 



Bn{z,a,b) :-- 



w{y)Pn{y, a, b)Pn-i{y, a, b) 



/i„_i(y - z) 



b 



1 c ^i^) - ^(y) 



+ 7 / Pniy,a,b)Pn-i{y,a,b)w{y)dy. 

K-i J a z-y 



A direct calculation produces two fundamental supplementary (compatibility) conditions 
valid for all z and these are stated in two different forms in the next lemma. 

Lemma 2. The functions y4„(2;,a, 6) and Bn{z,a,b) satisfy the conditions: 

B.n+i{z, a, b) + Bn{z, a, b) = {z - an)An{z, a, b) - 1/(2;), (^i) 
l + {z- an){Bn+i{z, a, b) - Bn{z, a, b)) = /3„+iA„+i(z, a, b) - /3„A„_i(z, a, b). {S2) 

It turns out that there is an equation which gives better insight into the a.„ and (in if {Si) 
and (5*2) are suitably combined. 
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Lemma 3. The functions An{z,a,b) , Bn{z,a,b) and 



n-l 



^Aj{z,a,b) 

j=0 



satisfy the condition: 



n-l 



B'^{z,a,b) + \/{z)Bn{z,a,b) + '^Aj{z,a,b) = (3nAn{z, a,b)An-i{z, a,b). (S*^). 

j=0 



3 Studies of Gaussian Unitary Ensembles 



It is the aim of this section to show that a quantity associated with Prob(n, a, b) defined in 
(12.31) satisfies a nonhnear pde for GUE via the ladder operators introduced in the previous 
section. To this end, we recall that in the present case w{x) = and / = M. Hence, an 
appeal to lemma 1 gives 

A„(^,a,6) = ^ + ^ + 2, (3.1) 
z — b z — a 

B^{z,a,b) = ^^ + ^, (3.2) 
z — b z — a 



where 



_ f,^ P^ib,a,b) 

hn{a,b) 
_ 2 P^{a,a,b) 
hn{a,b) 

_ _fe2 P„(&,a,&)P„_i(&,a,6) 
rn,b — — e - -, TT , [o.oj 

hn^i{a,b) 



_„2 Pn{a, a, b)Pn-i{a, a, b) 

y — Q ;^ 

hn-i{a,b) 

Substituting these into (5*1) implies 



(3.6) 



Rn,a + Rn,b — '^CUn, (3.7) 

r„+i,6 + rn,b = (b- an)Rn,b, (3.8) 

rn+l,a + rn,a = (a - an)Rn,a] (3.9) 
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while into {S2) , 



2 '^n,b'^n,a 

h — a 

J ^n,b^n,a 

a — b 



n-l 
j=0 

E 

j=0 





^ '^n,a '^n,b 

~ 2 2 2 ' 


n,a 


PnRn,aRn—l,a} 


'n,b 


PnRn,bRn—l,b} 


Rj,b 


= Pn 


Rn,bRn—l,a ~l~ Rn—l,bRn,a 


b — a 


Rj,a 


= Pn 


Rn,aRn—l,b ~l~ Rn—l,aRn,b 


a — b 



+ 2{R 

n-1,6 n,by 

Rn,a, 



The sum of the last two equations gives 

n-l 

n— 1,6 ~l~ Rn,b ~l~ Rn—l,a ^~ Rn,a) 

j=0 



,2 



2 

n,6 



^n,a ^n,b 



(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 



(3.15) 



where we have made use of (13. lip and (I3.12p to eliminate the terms Rn-i,a , Rn-i,b and 
arrive at the last equation. 



3.1 Toda evolution 



Before coming to the derivation of pde, we first establish some connections between /3„ , 
(n, a, b) and the auxiliary constants r„ ^ , r^^ , Rn,a and Rn^ , which will lead to a two 
variable analogue of the Toda equations for the recurrence coefficients q;„ and . These 
relations will be helpful in our further analysis. 

We start with taking the partial derivative of /i„ with respect to b in (12. 6p . which gives 

dbK = e-''Pl{b,a,b), (3.16) 

and consequently 



dh In hr 



-Rn.h- 



Since (3n = hn/hn-i , we find, 



1 

Jn 



dbPn — Rn-l,b — Rn,b, 



(3.17) 



(3.18) 
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and a similar computation yields 

■^da/3n = Rn-l,a — Rn,a- (3.19) 

Pn 

Again taking partial derivatives with respect to b on both sides of the equation 

= / Pn{x,a,b)Pn-i{x,a,b)e~'^ dx 



produces 

rb 

= Pn{b,a,b)Pn^i{b,a,b)e~^'' + / [9;,Pi(n, a, H ] P„_i(x, a, 6)e~'''rfx 

J a 

= Pn{b, a, b)Pn-i{b, a, 6)e"*'' + dbPi{n, a, b)hn-i{a, b), 
and consequently 



_yiPnib,a,b)Pn-i{b,a,b) 



9bPi(n, a, b) = -e 7 — , n ' = ^n,6- (3-20) 



A similar computation yields 



9„p,(„,a.(,) = e-- ^"'"-"-"'^"-;';-"-''' = r„,„. (3,21) 

ft„_i(a,6) 

Now, we are ready to prove: 
Lemma 4. For GUE, we have 

-^(5a + 56)/3n = 2K„i-aO, (3.22) 

Pn 

[da + di,)an = 2{(5n - Pn+l) ' 1- (3.23) 

Proof. Taking a sum of fl3.18p and 03.191) . the equation f l3.22p is immediate from fl3.7p . To 
estblish (13.231) . we note from the sum of fl3.20p and fl3.2ip that 

{da + 5b)Pl('^, a, b) = rn^a + = 2/3„ - n, 

where the last equality follows from fl3.10p . This, together with 02.110 . gives us (I3.23p . □ 
The equations 03.22p and 03.230 constitute a two variable version of the Toda equations. 
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3.2 Derivation of partial differential equation 



To this end, we set 

Hn{a,b) := ida + db)\nDnia,b), (3.24) 

where Dn is the Hankel determinant defined in fl2.5p associated with GUE. It is the aim of 
this section to derive a pde satisfied by if„ . Our strategy is to construct a system of partial 
differential equations in the functions r„_a , , Rn,b and Rn^a , because they provide a 
direct link to Hn ■ 

To see this, note that 

n-1 

In Dn{a, b) = '^lnhj{a,b) , (3.25) 

i=o 

and thus 

n— 1 n— 1 

da\nDr,{a,b) = ^da\nhj{a,b) = - ^Rj^a, (3.26) 

j=0 j=0 
n—l n—1 

dbln Dn{a,b) = J2dblnhj{a,b) = -J^Rj^b- (3.27) 

j=0 j=0 

Recall that we have from (13.71) 

Rn.b ~l~ Rn,a ^Qf^, 

and that from (I^TTD 

Q;„(a, b) = Pi{n, a, b) — Piiji + 1, a, b). 

Therefore we have that 

n— 1 n—l 

{Rj,b + Rj,a) = 2 ^ a,{a, b) = -2p^{n, a, b). (3.28) 

i=o i=o 

A combination of flX^ . flS:^ and flXTTD gives us 

n-1 

Hn = {da + db) \nDn = - J^(i?i,a + Rj,b) = 2pi(n, a, b) 

j=o 

= 2arn,a + 26r„,fe - 2^ - 2^ - 2/3„(i?„,, + R^m), (3.29) 
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where the second equahty of fl3.29p follows from fl3.15p . In view of (I3.20p and f l3.2ip , we note 
that 

daHn = 2rn,a, dbH^ = 2r„,b. (3.30) 
From fl3J9|) . f l318|) . (1311]) and fl3J2|) . we have 

daPn = f3nRn,a, (3.31) 

-fin, a 

dhf^n = PnRn,b- (3.32) 

This, together with f l3.10p . implies 

25bK,a + rn,b) = ~ (,2 + 2 / 



Expressing rn,a and rn,,a in terms of the partial derivatives of Hn , we have 

dlH,, + dadtHn = iMlil! _ (2n + daH^ + d,Hr,)Rr.,a, (3.33) 

dlH„, + 9,9bi/„ = IM^ _ (2n + daH^ + d,H^)Rn.b. (3.34) 

We may consider fl3.33p and fl3.34p as quadratic equations in Rn^a and Rn,b- Solving for 
them and substituting into (13. 29 p . we find, after some simplification, 

2hdbHn + 2adaHn - 2F„ =^{dlKn + dad^Hn)^ + A{daHnf{2n + daHn + 9,i/„) 

- ^(92/7, + a„abi/„)2 + 4(a6i/„)2(2n + daH^ + (3-35) 

(Note that in the above the signs of the square roots are determined by the signs of Rn^a 
and Rnfi ■ The former is positive and the latter is negative and the term 2n + daHn + di,Hn 
is also positive since it is the same as 4/3„ .) After clearing the square roots, we obtain that 
the function Hn defined in ( I3.24p satisfies the following pde: 

{{2hd,Hn + 2adaHn - 2Hnf - {{dlH^ + dad.Hnf + 4(9,i/„)'(2n + daH^ + 0^) 
- {{dlHn + dad^^Hnf + A{dbHn)\2n + daH^ + 9,if„)))' 
= A{{dlHn + dad^Hnf + A{daHnfi2n + + d,Hn)) 

X ((9,2/7„ + 9,9fe//„)2 + A{dkHnf{2n + + d^Hn)) . (3.36) 
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In our approach, the end points a and b are the "times" , although they play a distinct role 
from those in a two variable generalization of Painleve IV 



Suppose Hn is independent of a , the equation fl3.36p reduce to 

{dlHnf = A{hd^H:n - Hnf " ^{d,Hnf{2n + d,Hn), (3.37) 
which actually is the Okamoto-Jimbo-Miwa a -form of the Painleve IV equation [25j: 

{a"f = A{za' - af - A{a' + u^){a' + ui){a' + u^) (3.38) 

with 

uq = 2n, ui = 1^2 = 0. (3.39) 
The same conclusion holds if Hn is independent of b . 

Finally, it may be obvious, but worth pointing out that the solution if„ to the pde along 
with initial conditions does indeed yield the desired probability. This is because 

In Prob(n, a, 6) = /" Hn{t,t + b - a)dt + \nProh{n,0,b - a). (3.40) 
Jo 

3.3 Asymptotic independence of the extreme eigenvalues in GUE 

As an application of the pde derived in section 13. 2^ we show in this section that the extreme 
eigenvalues of GUE, when suitably centered and scaled, are asymptotic independent, i.e., 

lim Prob (^n, (^A^in + V2n^ n^^^/c > x, (y2n - A^ax) n^^^/c < 

= lim Prob (^(^A^^in + n^^Vc > lim Prob (^(^^2^ - A^^ax) n^^'^/c < , 

where Amin (Amax) denotes the smallest (largest) eigenvalue. 

Our method is to scale a and b near the edges of the spectrum and compute asymptotically 
the resulting pde. For this purpose, we let 

a = -v^ + c-^, 6 = v^-c-|^, (3.2) 



with c > , and note that 



^1/6 ^1/6 

da = d^, db = dy. 

c c 



Recall 

Hn{a, b) = (da + db) In Dn{a, b) = {da + db) In Prob(n, a, b), 



11 



and in the x , y variables this becomes 



^ Hn(-\^ + c^^,V2n-c-^] = {d^-dy)lnDn. (3.3) 



Let 



^(^' ^' •= + c;^' ^ - c ) • (3.4) 

After substituting the change of variables (13. 2p . the leading term of (I3.35P is of order n^^^ 
and produces the limiting pde 

- 8V2c^HdyHd,,H + 8V2c^y{dyHfd^H - 4{dyHfd^H + d^H{d'^H - d^dyHf 

+ dyH{8V2^x{d^Hf + A{d^Hf + {d^dyH - dlHf) = 0. (3.5) 

To ascertain whether the scaled smallest and largest eigenvalues may be described by their 
respective Tracy-Widom law for the extreme eigenvalues, we make use of a factorization 
ansatz, 

H{x,y) = f{x)+g{y), (3.6) 
where f{x) and g{y) satisfy the cr -form of a particular Painleve II. That is, 

\{nx)f = 2V2c'f{x)nx) - 2^c'x{nx)f - {nx))\ (3.?) 

\{g"{y)f = 2V2c'g{y)g'{y) - 2^c'y{g'{y)f + {g'{y)f. (3.8) 

We use this ansatz because if x = — oo or y = oo, then the functions / and g yield the 
correct solution of the pde and we believe that the solution should be their sum. Indeed this 
is true. An simple computation shows that (13.51) is satisfied identically. 

Now, we set 

P(x, y) = lim Prob (n, (^A^in + v^) n^^^c > x, (y2^ - A^^ax) n^^^c < y^ . (3.9) 
Recall that 

{d^ - dy)\nP{x,y) = lim H{x,y,n). (3.10) 

n— >oo 

Therefore the general solution of P{x,y) is of this form 

P(x, y) = F{x)G{y) exp [^{x + y)] , (3.11) 

where 

F{x) = lim Prob (^(^A^in + v^) n^^yc > , (3.12) 

G{y) = Jim Prob ((v^- A^ax) n'^Vc < , (3.13) 
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and \1/ is an arbitrary function. In view of fl3.ip . it remains to show \1/ = . To see 
this, note that 

hm F{x) = 1, hm F{x) = 0, (3.14) 

x—^—oo x—^oo 

hm G{y) = 0, hm G{y) = 1, (3.15) 

and 

hm P{x,y) = l. (3.16) 

x—>—oo,y—^oo 

Take a fixed z and let x + y = z . We see that 

hm P{x,z-x) = l = exp[^{z)], (3.17) 

for ah 2; . Hence \1/ = . 

An operator-theoretic proof of the asymptotic independence which also provide the rate of 
convergence to the factored Tracy- Widom distributions can be found in [10] . 



4 Studies of Laguerre Unitary Ensembles 



This section is devoted to the study of LUE. Hence, it is understood that all the notations 
hn{a,b) , an, /3n , Pi(n, a, 6) , etc. in this section are now associated with Pn{z,a,b) defined 
in (12.61) with w{x) = x"e~^ and I = (0, 00) . We will apply a similar theme to the LUE 
case as in the GUE case. 

By applying lemma [1] to the Laguerre weight w{x) = x°'e~^ , it is readily seen that 

. / ,N Rn Rn,a Rn,b f a i \ 

An{z,a,b) = \ \ -, (4.1) 

z z — a z — b 

Bn{z,a,b) = \ \ -, (4.2) 

z z — a z — b 
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where 



^« = irri^ I ' ^niy, &)y°"'e-^rf2/, (4.3) 

hn{a,b) 7„ 

= -'^-'^f^, (4.4) 

K,= -i-e-^^^, (4.6) 



a r 

^" = I / ^-(^' ^)Pn-iiy, a, b)y''-'e-ydy, (4.6) 

_ ^^^_^ Pn{a,a,b)Pn-i{a,a,b) 

hn-i[a,b) 
^ _,.^-> f„(M,6)P„-.(6,a,6) 
K-i[a,b) 



Substituting the above formulas into {Si) , we obtain 



From (^2) , we find, 



■,a ~l~ Rn,b 1) 


(4.9) 


n + Tn+i = a- a„-R„, 


(4.10) 




(4.11) 


, + rn+i,b = {b- an)Rn,b- 


(4.12) 




(4.13) 


^n,a ~ PnRn—l,aRn,ai 


(4.14) 


^n,6 ~ PnRn-\,bRn,b-i 


(4.15) 



n-l 



2 ■ 2 — + rn + a — - + a—^ + > Rj 

a b a b ^ 

j=0 

— R ( -^"-l."-^" Rn-lRn,a Rn-l,bRn + Rn-lRn,b\ , „n 

~a ^ b /' ^ ' 



n-l 

+ 2 + rn,a - a h > Rj,a 

a a-b a ^ 

3=0 

Rn—l,aRn ~l~ Rn—lRn,a . Rn—l,aRn,b ~l~ Rn—l,bRn,a 



Pnl - I - 7 I, I*--'- 'J 

\ a a — b / 
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2—- + 2-—- + r„,, - «^ + )^ R,,b 



j=0 

n f Rn-l,bRn + Rn-lRn,b Rn-l,aRn,b + Rn-1, bRn,a \ / /I 1 o\ 

= H b + b^a )■ ^^-^^^ 

From (S'2) , we find 

1 + r„+i - r„ + rn+i,a - rn,a + rn+i,b - rn,b = 0, (4.19) 

an{rn - rn+l) = /3„+li?„+l - (3nRn~l, (4.20) 
(a - a„)(r„+i_a - r„,a) = l3n+lRn+l,a - /3nRn~l,a, (4-21) 
(6 - an)irn+l,b - rn,b) = f3n+lRn+l,b " /3„-Rn-l,6- (4.22) 

The sum of (Hl6l) - (l4:T8|l yields 

rn + rn,a + Tn^b + ^ = 0, (4.23) 

where we have made use of fl4.9p . This equation can also be obtained by a telescopic sum of 
flCT]) . Summing flCTjl - flCTD . we see from gS]) and fl4:23|) that 

an = a + aRn,a + hRn,b + 2n + 1. (4.24) 
By fHTr|) - flTT5|) . it is easily seen that 

n— 1 n— 1 

i=o i=i 

= /3n ^(-Rri-l,a + Rn-l,b)Rn + (-Rn,a + Rn,b)Rn-l + Rn-l,aRn,b + Rn-l,bRn,a^ 

- 2r„(r„,a + r„,b) + (a - a)r„,a + (a - 6)r„,b - 2r„,ar„,6. (4.25) 

Now, we use (14.91) and (I4.23P to eliminate r„ and Rn ( Rn-i ), (14.141) and (14.151) to eliminate 
Rn-i,a and Rn-i,b in the above equation, it follows that 

n-l n-1 ^2 ^2 ^ 

O^-Rj.a + b''^Rj,b = (^n{Rn,a + Rn,b) + -^^(1 — i?„^fe) + "^^(1 " Rn,a) 
j=0 j=0 

+ {2n + a - a)rn,a + {2n + a - b)rn,b + 2rn,arn,b- (4.26) 



4.1 Toda evolution 

As in the GUE case, it is easily verified that, in the present case, we still have 

dalnhn = -Rn,a, dbln hn = -Rn,b- (4.27) 
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da In /3n — Rn-l,a " Rn,a, In — Rn-l,b " -Rn,fe, (4.28) 

with Rn a and -R„^f, defined in f l4.4p and fl4.5p . respectively, and 

da'Pi{n,a,h) = rn,a, dbPi{n,a,b) = rn,b- (4.29) 

daan = rn,a " ^^n+l.a, ^b^n = '^n.fe " '"n+l,?,, (4.30) 

where rn,a and r„^b are given in (14. 7p and (14.80 . respectively. 

With the above preparations, we are ready to state a lemma which gives a two variable 
version of the Toda equation for the recurrence coefficients a„ and Pn in the present case: 

Lemma 5. For LUE, we have 

{ada + bdb)l3n = /3n(a„_i - a„ + 2), (4.31) 

(ada + bdb - l)an = /3n- Pn+i. (4.32) 

Proof. From (I4.28p . it is easily seen that 

{ada + hdb)Pn ^ ^ _ ^^^^^^ ^ ^^^^^y ^4 33^ 

Pn 

Taking into account of (14.240 . this gives 

{ada + hdb)l3, 



Pn 

which is fOTD . 



- an + 2, (4.34) 



To establish ( I4.32p . we make the following decomposition 

{ada + &<9fe)a„ = ((a - a„)i9a + (& - a„)5;,)a„ + a„((9a + db)an. (4.35) 
By flCT]) and (g^S]), we have 

an(f^a + dh)(^n = an{rn,a + ?^n,6 " r„+i,a - r^+i^fc) = a„(l + r„+i - r^). (4.36) 

On the other hand, it follows from (001) . (|i2ID, (02ll and (S31) that 

((a - a„)9a + (fe - a„)<9b) 

= {a- an){rn,a - r„+i,a) + {h- an){rn,b - rn+i,b) 

= (3n{Rn-l,a + Rn-l,b) — Pn+l{Rn+l,a + Rn+l,b) 
= Pn ~ Pn+1 + Pn+lRn+1 ~ PnRn-l 

= Pn- Pn+1+ Oln{rn-rn+l), (4.37) 

where we also make use of (I4.20p in the last step of the above equation. Substituting (14.370 
and Km into gives us fICTD . □ 
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4.2 Derivation of partial differential equation 



We set 

Hn{a, h) := {ada + bdb) lnD„(a, b), (4.38) 

where is the Hankel determinant defined in (12.51) . It is the aim of this section to derive 
a pde satisfied by if„ . 

We note that, the equations (I3.26P and (13.271) still hold in the present case, it then follows 
from (K^ that 



71—1 n—1 



Hn = —a ^ Rj,a — b ^ Rj^b 
j=0 j=0 

= (a - « - 2n)r„^a + {b - a - 2n)r„^b - 2r„^„r„^b 



r. 



2 

- (3n{Rn,a + Rn,b) - ^(1 - i?,,,) - ^(1 - (4.39) 

In view of (I4.24p and (12. lip , we also have 

n— 1 n—1 

Hn = —a ^ Rj,a — b ^ Rj^b 
j=0 j=0 

n-1 

= ^(a + 2j + 1 - aj{a,b)) = n{a + n) + Pi(n,a,6). (4.40) 

j=0 

This, together with (I4.29p . implies 

daHn = rn,a, ObHn = r^fi- (4.41) 



Next, we derive representations of /3„ , Rn,a and Rn,b in terms of Hn and its partial deriva- 
tives. To this end, we use (I4.23P and (14. 9 p to eliminate r„ and Rn in (I4.13p . and then use 
(I4.14P and (I4.15P to eliminate the resulting Rn-i,a and Rn-i,b , it follows that 

— Pn{Rn,a + Rn,b) — "^^(^ ~ Rn,b) " ~ Rn,a) 

= -Pn + {n + rn,a + Tn^bf + a{n + rn,a + ^n^b) " T^a - T^b- (4.42) 

Inserting the above equation into (I4.39p . we obtain after some simplification that 

j3n = n{n + a)- Hn + arn,a + brn,b, (4.43) 
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or equivalently, taking into account of (14.411) . 

Pn = n'^ + an - Hn + adaHn + bdhH^ 
From f08|) . f H?Ti|) and f HlSj) . we further have 



n,b 



Using (I4.44P and ( I4.4ip in f l4.45p . it is readily seen that 



n +an- Hn + adaH^ + bdbHn)Rn,a 



n + an- Hn + adaHn + hdhHn)Rn,b- 



Solving the above quadratic equations for i?„ „ and Rn^b , we obtain 



Rr, 



R. 



-{adlHn + hdad,Hn)±V^l 

2{'n? + an — Hn + adaHn + bdbHn) ' 

-{bdiHn + adadM±VA-2 
2{n^ + an — Hn + adaHn + bd^Hn) ' 



where 



Ai = {ad'^Hn + bdad,HnY + A{daHny{n' + an-Hn + adaHn + bd^H, 



and 



(4.44) 



(4.45) 



(4.46) 
(4.47) 



(4.48) 
(4.49) 

(4.50) 



A2 = {bd^Hn + adad^HnY + 4(9,i/J^(n^ + an-Hn + adaHn + bd.Hn). (4.51) 



Finally, substituting (lOTj) . dOij) . flCTj) and (ICTj) into (ICTj) yields 

2(n^ + an - i:f„ + aS^i^n + bdbHn){Hn + {2n + a - a)daHn + {2n + a- b)dbHn + 2daHndbHn) 
+ {adlHn + bdadbHn){bdlHn + adadbHn) 

= 2{n^ + an-Hn + adaHn + bdbHn){±^i± v^) ± \/A^, (4.52) 

where Aj , i = 1,2 is given in f l4.50p and fl4.5ip . respectively. Denote by 

k := 2(n^ + an - if„ + adaHn + bdbHn){Hn + (2?t, + a - a)daHn + (2n + a - b)dbHn 
+ 2daHndbHn) + (a9^i/„ + bdadbHn){bdlHn + aS^Sfci/O 

and 

I := 2{n^ + an - Hn + adaHn + bdbHn), 
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we can rewrite equation (I4.52p in the following equivalent form: 

(((P - /2(Ai + A2) - - 4PA,A2{f + y^A^2 + Ai + A2))' 

- 16/6(AiA2)'(Ai + A2))' - 1024P(A^A2)^ = 0. (4.53) 



Suppose there is no a -dependence in if„ , the equation fl4.53p reduces to 

{bd^,H^y = {H^ + {2n + a~ h)d,Hr,f - A{d,H^f{n^ + an - + bdM- (4.54) 

The equation fl4.54p is nothing but the Okamoto-Jimbo-Miwa cr-form of the Painleve V 
equation [25] : 

(za"y = (a - za' + 2{a'f + (z/q + z/i + z/2 + 1/3)^^')' (4 55) 
- 4((t' + z/o)(a' + z/i)(fr' + z/2)(a' + z/3), 

with 

z/q = n, ui = n + a, 1^2 = J^s = 0. (4.56) 
We have the same conclusion if there is no 6 -independence in Hn ■ 



4.3 Scaling of PDE 

In this section, we will scale the pde obtained in ( I4.52p and show that its solution is asymp- 
totically equal to the sum of Tracy- Widom left and right distributions under certain delicate 
scaling. 

For this purpose, we set 

a = /3n, 13 >0, (4.57) 

and denote by 

L := 2 + /3-2v/2 + /3, R := 2 + (3 + 2^/2 + p. (4.58) 

Note that LR = . We then scale the left and right soft edges a and b as follows: 

a = Ln + cL'^^^n^^^x, b = Rn + cR^l^n^/^y, (4.59) 
with c > . Clearly, it is easily seen that 

dx dy 



and 

ad, 



•a = (^^'/' + ^) bd, = l^^n'/' + dy. (4.61) 
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Recall 

HrXa, h) = {ada + bdb) lnD„(a, b) = {ada + bdb) In Prob(n, a, b), (4.62) 
and in the new variables x , y this becomes 



Hn{Ln + cL^l^n^l^x, Rn + cR^'^n^'^y) 



(4.63) 



Let 

H{x, y, n) := -—HJLn + cL'^^^n^^^x, Rn + cR^'^n^'^y). (4.64) 

After substituting the change of variables f l4.59p . the leading term of fl4.52p is of order n^/'^ 
and produces the following limiting pde: 

4c Vl + PP^'^d^HdyH{H - xd^H - ydyH) 

+ dML"\d.dyHf + p"\dlHf + 2f3'/'L'/'d.dyHd'yH + Ac^^'' L''\dyHf) 

- dyHiR^'^d^dyHf + l3^/\dlHf + 2/3^/^R^/^d^dyHdlH + Acl3^'^R^/\d^Hf) = 0. 

(4.65) 



As in GUE case, it turns out that pde f l4.65p admits the following factorization ansatz: 

H{x,y) = f{x)+g{y), (4.66) 
where / and g satisfy the a -form of a particular Painleve II, respectively. More precisely, 

{nx)f = Ac'y^f{x)fix) - Ac'y^xifix)r - -^(/'(x))^ (4.6?) 

{g"iy)f = -Ac'y^giy)g\y) + Ac' ./TT^y{g\y)Y - -§j-,{9' {y)f ■ (4.68) 

It is readily seen that fl4.66l) - (l4.68l) satisfy (I4.65P identically. This proves that 

lim -^Hn{Ln + cL'^/'n^/'x, Rn + cR^'^n^'^y) = f{x) + g{y), (4.69) 

where / and g is Tracy- Widom left and right distribution, respectively. 
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